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Abstract 

Hot monomers are particles having a transient mobility (a ballistic flight) 
prior to being definitely absorbed on a surface. After arriving at a surface, 
the excess energy coming from the kinetic energy in the gas phase is dissi- 
pated through degrees of freedom parallel to the surface plane. In this paper 
we study the hot monomer-monomer adsorption-reaction process on a contin- 
uum (off-lattice) one-dimensional space by means of Monte Carlo simulations. 
The system exhibits a second-order irreversible phase transition between a re- 
active and saturated (absorbing) phases which belong to the directed perco- 
lation (DP) universality class. This result is interpreted by means of a coarse 
grained Langevin description which allows as to extend the DP conjecture to 

transitions occurring in continuous media. 
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I. INTRODUCTION 



Interacting particle systems are relevant to wide-ranging phenomena in physics, chem- 
istry, biophysics, ecology, etc. The concept of 'particles' is used in a broad sense, that is 
'particles' can be atoms, molecules, spins, individuals, etc, and whilst attention is drawn to 
the interactions among particles no attempt is made in order to achieve a detailed description 
(e. g. quantum mechanical) of the particle itself. Therefore, due to interactions, the oc- 
currence of complex behavior, such as phase transitions, self-organization, chaos, bistability, 
etc, may be observed0. 

Within this context, an active field of research is the study of far-from-equilibrium reac- 
tion systemsiH. Irreversible phase transitions (IPT) between active regimes (where reactions 
are sustained) and absorbing states (where reactions are no longer possible) have been re- 
ported in a great variety of models such as the Ziff, Gulari, and Barshad (ZGB) model 
for the catalytic oxidation of COi, the dimer-dimer mo del, the contact proces^, forest- 
fire mo delsi, etc (for a recent review see e. g.El). According to the Janssen-Grassberger 
conjectureili, irreversible reaction systems that exhibit a phase transition to a single absorb- 
ing state characterized by a scalar order parameter belong to the Directed Percolation (DP) 
universality class. This conjecture, stated a long time ago for unique absorbing states, has 
been further generalized for the cases where such states are non uniqueiS. A special case 
corresponds to non-equilibrium systems where, provided an IPT exists, there is in addition a 
local or global conservation of particles modulo two, such as the branching and annihilating 
random walks with an even number of offspringsS'lll. In these cases a new universality class 
emerges, commonly called parity conserving (PC) class, which is due to the existence of two 
statistically equivalent absorbing states at the critical pointlli. However, global conserva- 
tion of particles of modulo two may also lead to exponents in the PC class only when local 
spontaneous annihilation (IX — > 0) is highly inhibited. Then, at a coarse grained level, the 
relevant surviving processes are those conserving parity. In other words, parity conservation 
can be restored at a coarse grained level00. A nice example where global parity conser- 
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vation still leads to DP exponents is given by Inui et. a\S3. It is clear in this case that 
spontaneous annihilation must be taken into account. 

IPT are studied largely by means of Monte Carlo simulations and mean-field approaches. 
Recent developments of field-theoretic renormalization group techniques have provided a 
new theoretical framework where non-equilibrium phase transitions can be studied^. These 
techniques are able to identify the relevant processes in a given universality class although 
the quantitative predictions are still poor. 

So far, most of the simulations have been performed using discrete lattices where each 
particle fills a single site on the lattice and neighboring particles may react with a certain 
probability. In contrast, our knowledge on the behavior of irreversible reaction systems in 
continuous media is rather poor. In order to stress some dramatic differences that may arise 
for a reaction system when it is simulated off lattice, let us consider the simplest case of 
the B + B irreversible reaction which proceeds according to the Langmuir-Hinshelwood 
mechanism. 

B{g) + S^B{a) (1.1) 
B{g) + B{a) 
B{a) + B{a) ^0 + 2^ 

where g and a refer to the gas and adsorbed phases, respectively, while 5* represents a 
site on the surface. At first, we assume that S-species adsorbed on nearest neighbor sites 
react with unitary probability (P^ = !)• If we used a discrete lattice, reactions would be 
sustained indefinitely, i. e. the system could not irreversibly evolve into an absorbing state. 
However, considering a continuous media, the random adsorption of B particles of finite size 
a causes the formation of several interparticle gaps of size smaller than a. So, in the infinite 
time limit (t oo) the sample becomes imperfectly covered by -B-species separated by 
small interparticle gaps. Reaction is no longer possible and the system becomes irreversible 
trapped into an absorbing state (infinitely degenerated). The maximum jamming coverage 
attained in one dimension is Bj ~ 0.74759, which corresponds to the so called car parking 
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problem EE. 

In this paper we show that by introducing the adsorption of species with transient mobil- 
ity in a continuous one-dimensional medium, it is possible to open a window where reactions 
are sustained. However, by tuning the external parameter which controls the transient mo- 
bility of the particles it is possible to irreversible drive the system into an absorbing state. 

It should be mentioned that the study of reactions of atoms in the gas phase possessing 
thermal energy with adsorbed atomic CO species on metal and semiconductor surfaces is 
a topic of current interest. In contrast to thermally activated reactions among adsorbed 
species, i. e., the so called Langmuir-Hinshelwood mechanism, these kind of reactions take 
place under far-from-equilibrium conditions. Consequently, the determination of the un- 
derlying mechanism as well as the understanding of the dynamic behavior is challenging 
work. Within this context, very recently Kim et a@ have reported experimental studies 
of the reaction of hot iZ-atoms with adsorbed D-atoms (for further experimental works see 
references 

It should be noted that from the theoretical point of view a number of related models for 
random sequential adsorption with diffusion0'@ and desorptionS have also been proposed 
and studied (for a review see alsolil). However, interest in such studies is addressed to the 
asymptotic approach to the jammed state. In contrast, in this paper our interest is focused 
on the irreversible critical behaviour of a reactive system. 

So this work is devoted to the characterization of such IPT in continuous media and 
it is organized as follows: section 2 gives the description of the model and the simulation 
technique. In section 3 we discuss the results while conclusions and remarks are presented 
in section 4. 

II. THE MODEL AND THE MONTE CARLO SIMULATION METHOD 

In this paper, we study a ID off-lattice adsorption-reaction model in which particles of 
size a undergo a balhstic flight just after depositing on the substrate. The system evolves in 
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time under the following local rules, (i) A position x is randomly selected on the substrate. 
If the interval [x — a /2,x + a/2] is empty, then the adsorption trial is successful, otherwise it 
is rejected. So, it is clear that double occupancy of positions is forbidden, (ii) Right after a 
successful adsorption trial, a random direction is selected (left or right). Then, the particle 
undergoes a ballistic flight in the previously selected direction up to a distance R from the 
adsorption position x, provided that no other already deposited particle is found along the 
way. (iii) If during the flight one particle hits another previously adsorbed particle which is 
already at rest on the substrate, the following alternatives can occur: (1) the annihilation 
{B + B ^ 0) occurs with probability P^- Then, both particles react and leave the system. 
(2) Particles do not react (with probability (1 — Pr)), and the flying particle is frozen at the 
collision point. 

The ballistic flight mimics 'hot monomer' adsorption, allowing the incoming particle to 
transform its energy into degrees of freedom parallel to the substratum. The length of the 
flight R is finite in order to account for frictional dissipation. The model has two externally 
tunable parameters, namely R and P^. For P^ = one recovers the 'hot monomer' random 
sequential adsorption modefl while for P = and P^ = one has the Id car parking 
problem^. 

In order to simulate a continuous medium on a digital computer, one actually considers 
a discrete lattice. However, each site of size a is subdivided into 2^^ different adsorption 
positions. This high degree of discretization has provided excellent results when compared 
with the exact analytic solution of a related problemil. 

Preliminary results show that the system can undergo continuous IPT between a station- 
ary reactive state and an absorbing state without reactions when varying the parameters. 
This can easily be tested by considering the case P,. = 1 and P = (P > 1) which gives 
an absorbing (reactive) state, respectively. It should be pointed out that continuous IPT 
are dominated by fluctuations. Consequently, in a finite system and close to the critical 
point, the stationary state of the reactive phase can irreversibly evolve into the saturated 
state (absorbing state). Due to this circumstance, the precise determination of both critical 
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points and critical exponents is rather difficult. However, this shortcoming can be avoided 
by performing an epidemic analysis. For this purpose one starts, at t = 0, with a a config- 
uration close to one of the absorbing states. It is clear that different absorbing states will 
normally differ in the density of monomers. It should be pointed out that the dynamical crit- 
ical behavior of systems with infinitely many absorbing configurations is expected to depend 



to be independent of it. From the set of possible initial densities, a value pn is particularly 
important, namely the stationary density of inactive particles which results after the natural 
evolution of the system in the subcritical region has finished. The value p„ is relevant, since 
only for this value the natural dynamical critical behavior emerges. Preliminary simulation 
results show that p„ depends on the parameter P^, but their values have not been included 
in this work for the sake of space. The dependence of the critical behavior on the set of 
initial densities is the subject of an ongoing investigation. 

Consequently, the initial state for the epidemic analysis is generated by the evolution of 
the system very close to the critical point until poisoning is achieved. After generating this 
stationary poisoned state, we remove one or two particles from the middle of the system in 
order to create a small active area where adsorption is now possible. It should be noted that 
an empty area is considered to be active if it is longer than or equal to a. Then, the time 
evolution of the system is analyzed by measuring the following properties: (i) the average 
amount of active area at time t, A{t)\ (ii) the survival probability of the active area at time 
t, Ps{t)', and (iii) the average distance over which the active area have spread at time t, D{t). 
Finite size effects are absent because the system is taken large enough to avoid the presence 
of active area at the boundaries. For this purpose a sample of lO^'cr is enough. Averages are 
taken over 10^ to 10^ different samples. Near the critical point, the amount of active area 
is often very small. Then, we improve the efficiency of the algorithm by keeping a list of 
the positions where there is active area. Time is incremented by l/a{t), where a{t) is the 
amount of active area at time t. The time evolution of the active area is monitored up to 
t = 10^. At criticality, the following scaling behavior holds: 




However, static critical behavior appears 
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A{t) oc t\ 



(2.1) 



Psit)^t-', (2.2) 

and 

D(t)oct^/2 (2.3) 
where 6, t] and z are dynamic exponents. 

III. RESULTS AND DISCUSSION 

Preliminary simulations show that for = 1 it is possible to achieve a stationary reactive 
state in the large R limit (i. e. for R> 2) while in the opposite limit {R < 1.5) the system 
becomes irreversible saturated by -B-species. In order to obtain a quantitative description of 
the IPT we have performed epidemic studies around the critical edge. Figure 1 (a - c) shows 
log-log plots of A, Ps and D versus the time t, obtained for different parameter values. The 
three plots exhibit a power law behavior which is the signature of a critical state. Using 
smaller (greater) i?-values we observe slight upwards (downwards) deviations in three plots 
which indicate supercritical (subcritical) behavior (these results are not shown for the sake 
of clarity). Then, the critical exponents obtained by regressions are: 

= 0.308 d= 0.004 5 = 0.165 ± 0.003 2/2 = 0.625 ± 0.002 (3.1) 

These values are in excellent agreement with the exponents corresponding to the DP univer- 
sality class in 1 + 1 dimensionsi'i. Recently, extended series expansions calculations^ have 
provided very accurate values for the DP critical exponents, namely 

r] = 0.31368(4) 5 = 0.15947(3) z/2 = 0.63261(2) (3.2) 

Therefore we conclude that the studied adsorption-reaction model on a continuous medium 
belongs to the DP universality class like many other systems already studied on 



discrete lattices . It should be noticed that the present model has infinitely many absorbing 
states, so as in the case of the dimer-dimer mo del the DP conjecture holds for non-unique 
absorbing statesi at least as long as the absorbing states can be solely characterized by the 
vanishing of a single scalar order parameter. 

We have also studied the case of imperfect reaction, i. e. Pr < 1. Figure 2 shows a plot 
of the phase diagram. The phase boundary curve was determined by means of an epidemic 
analysis, as is shown in figure 1. The obtained critical exponents are 

r/ = 0.312 ±0.004 5 = 0.157 ± 0.003 2/2 = 0.631 ± 0.001 (3.3) 

Once again, these exponents are in good agreement with those corresponding to DP. Scan- 
ning the whole critical curve we obtain second order IPT's that belong to the DP universality 
class. However, the special case R ^ 00 merits further discussion. For Pr = 1 {Pr = 0) 
the system evolves towards a reactive (absorbing) state, respectively. Then, a transition is 
expected at some intermediate Pr value. In this case, the time evolution of the active area 
can be described by means of a mean field equation. In fact, the active area A{t) will grow 
(decrease) proportionally to A(t)Pr {A{t){l — Pr)), respectively; so 

^ = A{t)Pr - Am - Pr) (3.4) 

which leads to 

A(t) = Ao e(2^'-i)* (3.5) 

Therefore, = 1/2 is a critical value such as for Pr > 1/2 (P,. < 1/2) the active area will 
increase (decrease) exponentially in time, while just at criticality A{t) will remain constant 
{A{t) = Aq), which is consistent with a mean field exponent timf = 0. The predicted 
behavior is confirmed by means of simulations as it is shown in figure 3. By means of linear 
regressions the following exponents are obtained for P^ = ^: 

r/^0.0 S^l.O z/2^1.0 (3.6) 



8 



Then, our mean field estimate for t] is in good agreement with the simulation results. Re- 
grettably, we were unable to derive the mean field values for the remaining exponents. 

We conclude that the particular point P,- = 1/2, -R — > oo is a first order point (see figure 
2) which is not in the DP class which characterizes the whole critical curve. 

In the following, we give theoretical arguments by means of a coarse-grained Langevin 
description that support the result concerning the universality class of the model. First, 
note that the normalized variables needed to characterize the configurations of the system 
are the amount of active area a{x,t), the number of monomers in the system n{x,t) and 
the amount of inactive area v{x,t). These variables are not independent since we have the 
constrain a{x,t) + n{x,t) + v{x,t) = 1. It is clear from the above discussion that the time 
evolution of the system ends when a{x,t) = 0. Since a{x) ^ at criticality, this quantity 
can be chosen as the order parameter of the system. Then, we will try to describe the time 
evolution of the system near the critical point by means of two coupled Langevin equations, 
for instance, one for a{x,t) and the other for n{x,t). Due to the nature of the absorbing 
configurations, each term of these equations must vanish when a{x,t) 0. 
Let us consider the microscopic processes which are relevant to characterize the critical 
behavior of the system. First of all, both diffusion of a{x) and n{x) can be interpreted as 
successive adsorption- reaction processes. Within a one site description, both n{x) and a{x) 
will increase proportional to a{x). The reaction processes will contribute to the equations 
with a coupling term proportional to a{x)n{x). It is also clear that monomer flights will 
introduce terms proportional to a(x)^, a{x)^, etc. Since only the lower order terms are 
relevant for a renormalization group treatment0, we just keep the term proportional to 
a{xy. Then, we can write down the following Langevin equations: 

dn{x, t)/dt = ki'V'^a{x, t) + k2a{x, t) — k^a^x, t)^ — k^a^x, t)n{x, t) + rji{x, t) (3.7) 

9a(x, t)/dt = MiV^a(x, t) + U2a{x, t) — u^a^x, t)^ — u^ai^x, t)n{x, t) + 772(0;, t) (3.8) 

where r]i{x,t) and ■r]2{x,t) are uncorrected noises proportional to ^/a{x, t), ki and Ui are 
coefficients. This system of coupled Langevin equations is similar to that obtained for the 
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'pair contact processlaj which is one of the prototype systems with multiple absorbing states. 
Muiioz et all have shown that for large t the equation corresponding to the activity (equa- 
tion ( |3.8| ) for the present model) reduces to the Langevin representation of DP. Then, our 
simulation results are consistent with the above-presented theoretical arguments. The same 
authors have also shown that systems with many available absorbing configurations display 
strong memory effects that may lead to anomalous scaling. In addition to this, Mendes 
et aii have proposed a generalized hyperscaling relation which has proved to be valid in 
systems with multiple absorbing configurations. Simulation results on several lattice models 
with infinitely-many absorbing states support both theoretical argumentili@. The role that 
initial states play in the temporal evolution of the present model is under investigation. 



IV. CONCLUSIONS AND FINAL REMARKS 

A model for the irreversible adsorption-reaction of a single species on a continuous 
medium is studied by means of numerical simulations. We would like to stress and com- 
ment upon the following interesting features of the system: (i) in contrast to standard 
(reversible) transitions, non-equilibrium IPT can happen in one dimension, (ii) the stud- 
ied adsorption-reaction model clearly shows interesting new effects that may arise when a 
process is modelled on a continuous medium. Since the system always reaches a stationary 
reactive state when simulated on a discrete lattice but a final poisoned state can be observed 
on a continuous one (e. g. for Pr = 1 and R = 0), one may expect a crossover behaviour 
when the 'discretization degree' of the surface is tuned from one extreme to the other. This 
can be achieved by considering the adsorption on discrete lattices of species of arbitrary 
length r, i. e. r-mers. We found that the reactive random sequential adsorption (RRSA) 
of dimers (r = 2) always leads to a reactive steady state, whose stationary coverage is close 
to B ~ 0.5, and no poisoning is observed. However, the RRSA of trimers (r = 3) causes 
irreversible poisoning of the lattice with an average saturation coverage close to Qst ~ 0.7. 
In the case of dimers, two absorbing states of the type 
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...BBVBBVBBVBB.. 
.VBBVBBVBBVBBV.. 



(4.1) 



where V is an empty site, can be expected. So, during the RRSA the formation of several 
interfaces of the type ...BBVBBVVBBVBB... takes place. Due to coarsening we expect 
that in the asymptotic limit {t oo) both absorbing states will form two semi-infinite 
domains separated by an interface. The competition between these domains will keep the 
system in the reactive state for ever. Just by increasing r = 2 to r = 3, an infinite number 
of absorbing configurations appear in the system. So, the arguments developed above no 
longer hold and poisoning is observed. Consequently, an IPT is located between r = 2 and 
r = 3 and its precise location requires the study of the adsorption-reaction problem with 
particles of non-integer length. However, adsorption of r-mers of length r = 2 + e would 
cause the argument of the two competing absorbing states to fail. Therefore, we expect 
that the critical size, i. e. 'the discretization degree' is 2. (iii) To our best knowledge, 
this is the first off-lattice model which exhibits a second-order IPT in the DP universality 
class. Consequently, this results once again supports the DP conjecturei which can now be 
generalized to off-lattice models with infinitely many absorbing states. 

We expect that the interesting behavior of the present simple reaction model will stim- 
ulate further work on irreversible transitions and critical phenomena on continuous media, 
a field that, to our best knowledge, remains almost unexplored. 
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FIG. 1. a-c). Log-log plots of active area A{t), survival probability Ps{t) and epidemic diameter 
D{t) versus t. Results are obtained for a) Pr = 0.60 and R = 6.40, b) Pr = 0.80 and R = 1.685 
and c) Pr = I and R = 0.853. Averages are taken over (10^) runs with different poisoned initial 
states. 



FIG. 2. Phase diagram for 5 + B = reaction on the continuum (off-lattice). The active and 
adsorbing states zones are shown on a Pr versus R/ {R + 1) plot. The extreme values of R in the 
diagram are 0.853 and oo for = 1 and Pr = 0.5, respectively. According to that shown in Figs. 1 
and 3 the case given by the values R = oo and Pr = 0.5 is the only one which does not correspond 
to directed percolation universality class. 



FIG. 3. Log-log plot of the active area A{t), survival probability Ps{t) and epidemic diameter 
D{t) versus t for the singular case with R = oo. The solution of this case is exactly the one given 
by the mean field approximation. Dotted, full, and dashed lines correspond to epidemic studies for 
Pr < 0.5, Pr = 0.5 and Pr > 0.5, respectively. 
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